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(54) Signal processing system 

(57) A CMOS integrated signal processing system 
for a sanpling receiver includes a tinning recovery dr- 
cuit. wherein an on-chip numerically controlled oscillator 
is operative at periods T that are initially equal to the 
nominal baud rate of the signals controls a sine interpo- 
lator receiving samples at the sampling rate. A loop filter 
is coupled to the sine interpolator and to the numerically 
controlled oscillator. The arrangement is capat)le of 
handling various symbol rates. The system includes a 



circuit for carrier recovery, having a second on-chip 
numerically controlled oscillator, a digital derotetion cir- 
cuit responsive to the second numerically, controlled 
oscillator, accepting an in phase component and a 
quadrature component of the sampled signals. An 
adaptive phase erra estimation circuit is coupled in a 
feedback loop. 
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Description 



This invention relates to processing signals received from a communicatiorrs channel. More particularly this inven- 
tion relates to an integrated signal processing system for demodulating signals suitable for use in the transmission of 
television signals. 

Encoded transmission of inherently analog sigrials is increasingly practiced today as a result of advances in signal 
processing techniques that have increased the bit rate achievable in a channel. At the same time new data compression 
techniques have tended to reduce the bandwidth required to acceptably represent analog information. 

Various modulation techniques have been employed in digital communications. For example quadrature amplitude 
modulation (QAM) is a relatively sophisticated technique favored by practitioners of digital radio communications. This 
method involves two separate symbol streams, each stream modulating one of two carriers in quadratic e. A transmitted 
QAM signal can be represented by the equation 



is a finite sequence of transmitted symbols; 
g(t) is a real valued transmit filter; and 
T is the synrtool period. 

As will be apparent to those skilled in the art this is equivalent to the modulation of two real-valued baseband pulse 
amplitude modulated (RAM) signals by the carrier signals cos {a^X) and sin respectively. As used heran the first 
term in the above equation is referred to as the "in-phase" component, and the second term is called the "quadrature** 
component. 

This system achieves spectral efficiencies between 5 - 7 bits/sec-Hz in multilevel formats such as 64- and 256- 
QAM. QAM Is particularly useful in applications having a high signal-to-noise ratio. However double sideband nxxJula- 
tion is required, which requires increased channel bandwidth for the same symbol rate over single or vestigial modula- 
tion schemes. Furthernme cross-coupled channel equalizers are generally needed to cancel linear distortion in the 
channel, which adds to the overall complexity of the system. 

A variant of QAM is quadrature phase shift keying (QPSK), in which a signal constellation consisting of four sym- 
bols is transmitted, each having a different phase and a constant amplitude. The scheme is inplemented as the sum of 
ortiiogonat components, represented by the equation. 



where 0^ can be any of {0. 7c/2. n. 3jd2]. It is necessary to transmit both sidebands in ader to preserve the quadrature 
infbmriation. The QPSK modulation scheme has been adopted by the ITU-T as an International standard for direct dig- 
ital satellite broadcasting. In Europe 16-QAM and 64-QAM are used in the digital video broadcasting (DVB) standard 
for digital cable broadcasting. Both QAM and QPSK have similar coding schemes, which are generally described in Rg. 
1 with reference to MPEG transport layer packets, wherein QPSK and QAM are implemented according to the stand- 
ards DVB-S (European Telecommunication Standard PrETS 300 421) and DVB-C (European Telecommunication 
Standard PrETS 300 429) standards. MPEG is a standard well known to the art, In which data is g'ouped in a plurality 
of packets, each of which contains 188 bytes. This number was chosen for conpatibility with asynchronous transfer 
mode (ATM) transmissions, another known telecommunication standard. Various aspects of the coding process are 
specified in the respective DVB standards, including: randomization and sync inversfon for synchronization; Reed Solo- 
mon encoding; Forney interleaving; convolution encoding in the case of DVB-S, and byte to m-tuple mapping, and dif- 
ferential mapping in the case of DVB-C. 

The art is presentiy striving to more efficientiy transmit video and audio data in applications such as cable and 
direct satellite television using digitaf techniques. 

For a better understanding of these and other objects of the present invention, reference is made to the detailed 
description of the inventfon, by way of example, which is to be read in conjunctfon with the following drawings, wherein: 

Fig. 1 is a block diagram illustrating QAM and QPSK coding and modulation; 
Fig. 2 is a block diagram showing modulation and demodulation in a communfoation system; 
Fig. 3 is a more detailed block diagram of the tuner and I.Q denKxIulator shown in Fig. 2; 
Fig. 4 is a diagram illustrating a portion of tiie modulator circuit shown in Fig. 2; 



x{t)=V2 cos(w,t) Re{aJg(t-mT) -V2 sin(w,t) X lm{a„)g(t-mT) 




wherein 
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Fig. 5 is a plot shewing interpolated samples disposed between known samples: 

Rg. 6 is a block diagram of a carrier recovery and a timing recovery circuit in a receiver according to the prior art: 
Fig. 7 is a more detailed block diagram of the carrier recovery circuit shown in Fig. 6; 
Fig. 8 Is a more detailed block diagram of the timing recovery circuit shown in Rg. 6; 
5 Fig. 9 is a diagram that illustrates the operation of the Gardner algorithm: 

— ^ Rg. 1 0 Is a block diagram of a canier recovery and a timing recovery circuit according to the Inventkxi; 

Rg. 11 is a tvore detailed, partially schematic block diagram wtvch illustrates the tinung recovery circuit shown in 
Rg. 10: 

Rg. 12 is a schematic illustrating a numerically controlled oscillator in the tinrmg recovery circuit shown m Rg. 1 1 ; 
10 Rg. 13 is a diagram illustrating the operation of the circuit of Fig. 1 1 ; 

Rgs. 1 4a and 1 4b are plots of sine pulses having different delay values which are he^l in understancBng the oper- 
ation of the circuits illustrated in Rgs. 11-13: 

Rgs. 15a, 15b, 15c are plots which indicate the process of sine interpolation: 
Rg. 16 is a schematic of a filter used in the sine interpolator unit of the circuit shown in Rg. 1 0: 
15 Rg. 1 7 is a schenriatic of a linear interpolation unit that can be used \n the drcuit of Rg. 10: 
— ^ Rg. 1 8 is a schematic of a matched filter that is used in the circuit of Fig. 10; 
Rg. 19 is a block diagram of a canier recover circuit according to the inventbn: 

Rg. 20 is a k>lock diagram of an adaptive phase tracking circuit that is used in the carrier recovery circuit depicted 
In Rg. 19: 

20 Rg. 21 is an electrical schematic of the hopping adder used in the circuit illustrated in Fig. 19; 

Rg. 22 Is a tiock diagram illustrating a carrier recovery and a timing recovery circuit according to an attemate 

ennbodiment of the Invention: 

Fig. 23 is a more detailed diagram of a Hilbert filter used in the drcuit shown in Rg. 22; 
Fig. 24 is a diagram illustrating constellation rotatk>n error; and - 
25 Rg. 25 is a bk)ck diagram illustrating a second altemate embodiment of the inventk>n. 

It is a primary object of the present Inventkm to provkJe an improved signal processing system for the communtoa- 
tion of data in a constrained channel. 

It is another object of the invention to provide an improved, economical apparatus for receiving and decoding data 
30 at high bit rates, such as vkJeo and audio signals. 

It is a further object of the invention to provkie improved apparatus that economically and reliably provide locking 
for the demodulation frequency according to the modulator frequency in a communication system. 

It is still anotiier object of the invention to provide improved apparatus for economically and reliably locking tiie data 
sampling frequency according to the rate of transmitted data in a communication system. 
35 These and other objects of the present invention are attained by a signal processing apparatus for processing sig- 
nals that are sampled by a sampler operative at a sampling rate. The apparatus comprises a clock operative at the sam- 
pling rate, a first numerically controlled oscillator operative at periods T that are Initially equal to the nominal baud rate 
of the signals, an interpolator, preferably a sine interpolator, receiving samples at the sampling rate, and a loop filter 
coupled to the sine interpolator. The loop fitter has an output responsive to a difference between the periods T and a 
40 received symbol rate of the sampled signals. The first numerically controlled oscillator is responsive to the toop f ater and 
generates an output signal that is representative of an interpolatton distance between succeeding samples. The sine 
Interpolator interpolates the received samples according to the interpolation distance, and produces an output signal 
representative of the interpolated samples. 

In an aspect of the invention the input signal Is modulated, and the apparatus further comprises an l,Q demodula- 
45 tor. First and second analog-to-digital converters are respectively coupled to an in phase output and a quadrature out- 
put of the demodulator, wherein the sine interpolator accepts in phase and quadrature signals. 

In aripther aspect of the invention the loop filter accepts an in phase component of the interpolated samples and 
the error signal is computed according to the equation 

error(r) = l[r.|][l(rHI(r-T)J 

wherein I is the in phase component, T is the symbol period, and r is tiie interval between alternate samples. 
55 In yet another aspect of thte invention tiie loop fitter accepts an in phase conponent and a quadrature component 
of the Interpolated sanples tiie error signal is computed according to the equation 
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efror(r) = l[r-|]ll(r)+l(r-T)l + Q[r.2[Q{r)+Q{r-T)l 



wherein I is the in phase component, Q Is the quadrature component. T is the symtx)i period, and r ts the intaval 
between alternate sanples. 

In still another aspect of the invention the first numerically controlled oscillator, the sine interpolator, and the loop 
filter comprise an integrated semiconductor circuit, preferably a CMOS circuit. 

The apparatus includes a matched f ilter which has an input coupled to the sine interpol^or and an output coupled 
to the loop filter. Preferably the matched filter is a square-root raised cosine f Pter. 

The output of the first numerically controBed oscillator comprises a f rst output signal that is generated whenever 
the state n exceeds a division of the symbol period, and the sine interpolator gen^-ates an output in response to the 
first output signal. The output of the first niffnericaily controlled oscillator includes a second output signal that is repre- 
sentative of a value A accprding to the equation 



wherein system dock is the clock rate, baud rate is the nominal baud rate, and a is a state representative of a nun^er 
of elapsed operative periods of the first numerically controlled oscillator, and the sine interpolator emits an interpolated 
sample upon receiving the second output signal. 

Accading to an aspect of the invention the sine interpolator is a unit comprising a first sine interpolator that receives 
an in phase component of the sanrples. and a second sine interpolator that receives a quadrature component of the 
samples. Preferably the sine interpolator comprises a finite impulse response filter having a bank of coefficients. 

In yet another aspect of the invention the bank of coefffcients conprises a plurality of banks, and sine interpolator 
Is provided with an addressable memory containing a plurality of coeffidenta 

In another aspect of the invention the sine Interpolator performs a plurality of sine interpolations that precede and 
follow a required sine interpolation point, and the apparatus further comprises a linear interpolator that performs linear 
interpolation on the plurality of sine interpolations. 

The invention provkles a circuit for processing modulated signals, comprising a semiconductor Bitegrated earner 
recovery circuit operative to control a demodulator that includes a second numerically controlled oscillator, and a digital 
derotation circuit which is responsive to the second numerically controlled oscillator and accepts an in phase con^ 
nent and a quadrature component of sanrpled signals. The carrier recovery circuit further comprises a phase en-or esti- 
mation circuit coupled to an ou^ut of the derotation circuit and a loop filter coupled to an output of the phase error 
estimation circuit wherein the second numerically controlled oscillator is responsh/e to the loop filter. 

In a further aspect of the invention the circuit adaptively estimates the phase en^or and executes a least-mean- 
square algorithm. It comprises first and second slicers which accept a derotated in-phase value and a derotated quad- 
rature value respectively: first and second subtracters for respectively determining first and secorxJ differences between 
the derotated in phase value and the sliced in phase value, and between the derotated quadrature value and the slk»d 
quadrature value: and an angulator, accepting the first and second differences and outputting a phase error estimate. 

Preferably the integrated circuit is a CMOS circuit 

The invention provides a signal processing apparatus for processing modulated signals at a modulatfon carrier fre- 
quency, comprising a demodulator, a sampler operative at a sampling rate on an output of the demodulator, and a car- 
rier recovery circuit operative to control the demodulator in accordance with the rrexJulatfon canrier frequency. The 
earner recovery circuit comprises a second numerically controlled oscillator: a phase error estimation dreuit: a digital 
derotation dreuit responsive to the second numerically controlled osdilator and accepting an in phase component and 
a quadrature conponent of sampled signals; and a loop filter coupled to an output of the derotatfon circuit; wherein the 
second numerically controlled oscillator is responsive to the loop filter. The sairpier and the carrier recovery circuit are 
integrated In a semiconductor integrated dreuit preferably a CMOS dreuit 

In an aspect of the invention the carrier recovery circuit adaptively estimates the pfiase error according to the least- 
mean-square algorithm. It comprises f rst and second slicers which accept a derotated in-phase value and a derotated 
quadrature value respectively; first and second subtracters, for respectively determining first and second differences 
between the derotated In phase value and the sliced in phase value, and between thb derotated quadrature value and 
the sliced quadrature value; and an angulator. accepting the first and second differences and outputting a phase error 
estimate. 

The invention provides a signal processing apparatus for processing modulated signals that are sanpled by a sam- 
pler operative at a sampling rate. The signals have a nominal baud rate. The apparatus comprising a dock, operative 
at the sampling rate; a first numerically controlled osdilator operative at periods T that are initially equal to the nominal 
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baud rate: and a sine interpolator receiving samples at the sampling rate. A first loop filter, ooupled to the sine interpo- 
lator and the first nunnencally controlled oscillator, has an output responsive to a difference between the periods T and 
a transmitted symbol rate of the sampled signals, wherein the first numerically controlled oscillator is responsive to the 
first loop filter and generates an output signal that is representative of an interpolation cfstance between succeeding 

5 samples, and the sine interpolator interpolates the received samples according to the interpolation distance, ar^ pro- 
duces an output signal representative of the interpolated sanies. A carrier recovery circiit comprises a second numer- 
ically controlled oscillator, a phase error estimation circuit, a digital derotation circuit responsive to the second 
numerically controQed oscillator and accepting an in phase compcnent and a quadrature component of sampled sig- 
nal, and a second loop filter coupled to an output of the derotation circuit wherein the second numerically controlled 

w osciOator is responsive to the second loop filter. The first and second numerically controlled oscillators, the sine inter- 
poteita. the f rst and second loop titers, and the digital derotation circuit are integrated in a sennoonduetor integrated 
circuit 

In one form of the invention an analog-to-dlgital converter unit samples an input signal at the sampling rate and has 
an output coupled to the sine interpolata. The analog to digital converter unit is integrated in the semicondirctor inte- 
rs grated drcuit 

In still another aspect of the invention the input signal is modulated, and the apparatus further comprises an I.Q 
demodulator, and a sampler comprising first and second analog-tOKiigital converters respectiv^y coupled to an in 
phase output and a quadrature output of the demodulator, wherein the sine interpolator accepts in phase and quadra- 
ture outputs of the sampler. 

20 The organizatton of a digital receiver 50 is shown in Rgs. 2 and 3. Although the invention is explained with refer- 
ence to a particular digital receiver, rt will be understood by those skilled in the art that it can t>e practiced in many forms 
of modulation and demodulation wherein it is needed to accurately lock the demodulation frequency to the modulation 
frequency. Similarly the teachings herein are applicable to many systems wherein carefully synchronized conversion of 
an analog signal to a digital signal is required. 

25 Rg. 2 illustrates a block diagram of a eonvnunications system which includes a digital receiver 50. A modulator 20 
modulates arid transmits a signal across a communications channel 22, which is initially accepted by a receiva- front 
end 23, and demodulated in an l.Q demodulator 21. The demodulated signal is sampled in an analog-to-digital con- 
verter ADC 60. Timing recovery is performed in timing recovery circuitry 62. Carrier recovery is accomplished in a car- 
rier recovery circuit 64. The receiver typically includes an automatic gain control (AGC) circuit 66. The sampled data is 

30 processed through a matched filter 68, and is then subjected to a sllcer 69 and erra con-ection circuitry 72, as dis- 
cussed below. 

In the receiver front end 23. a radio frequency amplifier 52 is coupled to the comnujnications channel 22. The chan- 
nel is typically a constrained channel, such as a satellite downlink, or a television cable, but can be any comnujnications 
channel. The output of the radio frequency amplifier 52 is demodulated to a fvst intennediate frequency by the first 

35 demodulator 24. and then passed through a band pass filter 25. The first denxxiulator 24 is of a known type. Its fre- 
quency is controlled by a voltage-controlled osdilator 33, programmed through a microprocessor interface 29 operating 
through a digitel-to-analog converter DAC 37. The intermediate frequency (IF) signal output from the amplifier 31 is 
demodulated in the I.Q demodulator 21 (Fig. 2) to a complex baseband representation using demodulator subunits 1 , 
2 included therein, and passed through low pass fitters 3. 4 to produce the output s^nals kJata and Qdata. which are 

40 the in phase and quadrature components respectively These signals are converted into a digKal representation using 
a high speed analog-to-digital converter unit ADC 60. 

The signal received from the channel 22 has been modulated as shown in Fig. 4. Pulses comprising the in phase 
component Idata and the quadrature component Qdata are filtered through raised root cosine filters 5 and 6 respec- 
tively, modulated on orthogonal carriers at the can-ier frequency co^ in multipliers 7. 8, and summed in an adder 9. Typ- 

45 ically signals in adjacent passbands are transmitted ak>ng with the desired signal 10. In typical applications the excess 
bandwidth of the shaped pulses exceeds the Nyquist minimum by 35%. Optionally, ttie modulator may include an inter- 
mediate frequency (IF) stage (not shown). For a symbol period T of 33 ns. the Nyquist frequency is approximately 15 
MHz. equivalent to approximately 30 MHz in ttie passband. For such a signal, a 40.5 MHz channel at 3 db rolloff would 
be required. The adjacent channels, channel distortion, and noise all nust be conskfered in tiie design of the receiver 

so and demodulator. The details of demodulation depend on the communications channel. 

Referring again to Rg. 2. a high speed analog to digital converter. ADC 60. provides an output which is used for 
timing recovery circuitry 62, which insures accurate sampling by ADC 60. In order to successfully recover the original 
data, the receiver is required to aceonplish at least tiie foltowing tasks: 

55 (a) lock tfie demodulation frequency to the modulation frequency; 

(b) lock the data sampling frequency to the transmitted symbol rate; 

(c) adjust tiie gain of the tuner for optimal signal to noise; 

(d) complete the pulse shaping to minimize intersynrtool interference (ISI) 

(e) perform a Nyquist filter operation to reject out-of-channel noise; 
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A carrier recovery circuit 64 controls the I.Q denrxxlulator 21 such that the correct frequency and phase are recov- 
ered. The AGO circuit 66 feeds t>acK to radio frequency amplifier 52. Both the AGC circuit 66 and the carrier recovery 
circuitry 64 are coi^led to the output of ADC 60. The main digital data stream from the ADC 60 is f itered by a matched 
filter 68. whicti precisely matches the characteristics of the transmit filters (not shown). A sllcer 69 extracts the data from 
the filter output 68. determines the nearest legal constellation point, and applies a representation in appropriate format 
to the error correcting circuitry 7Z The specification for the error correction is specified in the DVB specification, and will 
not be further discussed as it is outside the scope of the invention. 

Before proceec£ng to a further detailed description of the preferred emtxxliment. it is believed tf^t the invention can 
be more clearly understood by comparison to Figs. 6 - 9, which illustrate a conventional approach to the prot^ems of 
carrier recovery and tirrang recovery. An off-chip l,Q quadrature demodulator 10 is employed in the tuner. An external 
voltage controlled oscillator 42 is oontroHed by carrier recovery circuitry block 44 in ord^ to lock the demodi^tor dock 
to the frequency and phase of the transmitted signal. Typically the voltage controlled oscillator 42 can be adjusted 
across a range of several MHz. The sampling of the transmitted circuit is locked to the transmitted symbol rate and 
phase by an external voltage controlled oscillator 46, which is typically a crystal oscillator. The oscillator 46 is controlled 
by timing recovery drcuitry 48. but is generally limited to a range of a few hundred KHz. Because different satellite trans- 
ponders and cable television systems use (Afferent symbol rates, a plurality of voltage controlled osdilators (not shown) 
may be required to accommodate the diversity of transmission arrangements, or the device imy be limited in its appli- 
cation to a particuiar transmission system. The demodulated data is sampled at T/2. the system clock rate in analog- 
to-digital converters 45. 47. At a symbol rate of 30 Mbaud. the system clock is therefore running at 60MHz. T/2 sampling 
is required for a Gardner timing recovery loop. However other drcuitry outside the timing recovery kx)p is generally 
ckxtod a the stower rate of T to simplify implementation and reduce circuit area. After filtering in matched f aters 54. 56 
the I.Q outputs 58, 59 are sent to a slicer arxl error correcting drcuitry (not shown) as required by a particular applksa- 
tion. 

The matched filters 54. 56 are typically implemented as square-root raised cosine matched filters, having an 
excess bandwidth a = 0.35 to conform to the DVB specification. These filters match a transmit f iter (not shown) that 
was errpbyed in transmitting the input signal so as to restore the signal to its pretransmisskxi character. 

The carrier recovery circuitry 44 may be implemented as a conventional Costas t.oop. referenced generally at 61 
in Rg. 7. After a channel change there may be a significant frequency error which has to be determined before phase 
can be acquired. A frequency-lock-loop of the type where the error signal 67 is proporttonal to the frequency error, or a 
frequency sweeping scheme can be used to determine the initial frequency en-or. 

The timing recovery circuitry 48 is shown in more detail in Fig. 8. A conventional Gardner algorithm is employed, 
which acquires the timing sanple point for Ty2 sampling by the A/D converters 45. 47, and the units 55, 57. Even-num- 
bered samples are used as data samples, while odd-numbered samples correspond to zero aossings. An error is com- 
puted according to the equation 



lis the in phase output; 

Q is the quadrature output; 

T is the symkx)! period; 

r is the sample time of the even sample. 
The Gardner algorithm is explained in greater detail in A BPSK/QPSK Timing-Error Detector for Sampled 
Receivers, Gardner, Floyd M.. IEEE Trans. Comms. COM-34, May 1986, pp. 423-9. Qualitatively, the error signal indi- 
cates the direction the timing sanple point has to move in order to lock odd-numbered T/2 samples to the midpdnt 
between samples, and to kxk even-numbered sample points to the near optimal sample point As the data is rand- 
omized in practical systems, there are an adequate number of zero crossings. 

Accading to the algorithm the difference is calculated between a succeeding even-numbered sanple and a pre- 
ceding even-numbered sample in subtracters 38. 39. This difference is multiplied in multipliers 41 . 43 by the odd-num- 
bered sample therebetween. The case where the sample point is too early will be explained with reference to Fig. 9. In 
the case of a failing edge 70. the intermediate odd-numbered point 71 has a positive value. The difference between the 
even-numbered points 73. 74 is negative, because the succeeding pdnt 73 is closer to the zero line than the preceding 
point 74. Thus the product has a negative value. For a rising edge 80, the odd-numbered intermediate point 78 has a 
negative value. The difference t>etween the succeeding even-numbered points 82, 84 is positive. Thus the product also 
has a negative value. 

A similar analysis reveals that in the case of late sampling, the Gardner algorithm yiekis a positive value in both 
rising and falling edges. In the interest of brevity it will not be repeated. 




where 
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The error value indicates in wTtich direction the sample point must move to align it correctly. 

When the sample point is correct the error value reported by the Gardner algorithm is zero, except for noise and 
Inter-symbd interterence (IS!) effects. However the noise effects have a zero mean. 

After processing, the even-numbered samples are sliced to give the reconstructed data which is applied to the error 
6 correcting circuitry. The odd-nurhbered samples are discarded before derotation. Of course it is equally possible to slice 
the odd-number sarrples amd discard the even-nurrtsered samples with appropriate nxxfif ications of the timing recovery 
circuitry. Similarty. sampling at other divisions of T coukj be employed. Although locking with odd-numbered samples at 
the zero crossing points does not guarantee that the optimum sampling point has been found for the data, the scheme 
works well in practice. Alternative schemes which more dosely lock onto the point having a maximum likelihood of oor- 
10 redness could be used, but these are generally more corrplicated to implement 

The logic in the timing recovery drcuit 46 operates at T because an error estimate is generated only every other 
sample. However it requires two T/2 cycles for an even-numk)ered sample to propagate to the next even-numbered posi- 
tion in the shift registers 63, 65, 75. 77 (Fig. 8). Thus the samples are docked through at T/2. 

A preferred embodiment of the inventk)n will now be explained initially with reference to Rg. 10. It has an advantage 
IS over the conventional circuit discussed above, in that the carrier and timing loops have been inrplemented entirely within 
the digital domain, and are integrated on a CMOS chip. As a result externa) system component costs are reduced. A 
further advantage of the anangement is that the denxxiulator can work at many different synfix>l rates, or with variable 
symbol rate technology As discussed above, the prior art solution required the variable crystal controlled osdilator to 
match the symbol rate. 

20 A conventional off-chip l.Q baseband demodulator 140 is employed. A suitak)ie I.Q quadrature baseband demodu- 
lator is the GEO Plessey SL1710 l.Q denxxiulata. The external sample timing recovery k>op has been replaced by a 
fixed frequency system clock 120. which docks the logic. It must at least equal the data Nyquist frequency or otherwise 
ensure that the data Nyquist frequency is met An on-chi p interpolator unit 130. which, under control of the timing recov- 
ery loop 125, gene rates synchronous T/2 spaced sample valu e. At each system dock tick, either one or zero T/2 sam- 
ps pies are generated by the interpdator unit 130. In the event that a sample has been generated, subsequent hardware 
nxxJules are so informed by the assertion of a Valid' contrd strobe 1 70 (Rg. 1 1). 

The external carrier recovery loop illustrated in Rg. 6, induding the external voltage controlled oscillator 46, has 
been replaced with a fixed frequency external crystal osdilator 145. The osdilator 145 cooperates with an on-chip dig- 
ital derotater 150, and an on-chip carrier recovery loop 155. Both the derotater 150 and the earner recovery loop 155 
30 are operative virith T-spaced samples. The I and Q outputs 1 52, 1 54 are applied to the dicer and error correction drcuits 
in accordance with the DVB specif icat'on. 

Timing Recovery 

35 The timing recovery circuit according to the invention is shown in greater detail in Rgs. 11 and 12. Quadrature 
demodulated data is sanpled at the system dock rate, which, as explained above, must be at least equal to Nyquist 
frequency of the input data. As illustrated in Rg. 1 2. an on-chip numerically controlled osdilator. shown generally at 21 0. 
keeps count of symbol time. The state n of the numerically controlled osdilator 210 represents a fixed point count of 
the nunrber of symbol periods which have elapsed. At eadi system dock tick, the state n increments by a value in a 

40 register 260 equal to the (nominal baud rate / system dock rate), adjusted from the nominal using a control signal 262. 
The reciprocal of the value in the register 260 is taken in a reciprocal generator 240. The reciprocal is m Jtiplied by the 
value [(20) mod 1 1/2 in a multiplier 241 , and limited to a value that is. less than 1 in saturation block 242. 

Referring now to Rgs. 1 1 and 13. an example is shown that illustrates the operation of the cffcuit shown in Rg. 1 1 . 
wherein the (nominal baud rate / system dock rate) = 0.4. Whenever the state O passes a T/2 marie, indicated by the 

45 \CM& row of upward directed arrows in Rg. 13, the numerically controlled osciDator 210 outputs a signal 215. and a 
value A 217 which represents an interpolation distance between successive sanples given by 



The signals 215. and A 21 7 are accepted by a sine interpolator unit 222. conprising an irxJividual sine interpolator 221 a. 
221b for each of the in-phase and quadrature components,. The interpolator unit 222 then generates a sanple value 
based on the interpdation distance. A has a value greater than or equal to 0 and less than 1. but is represented as a 
55 fixed point number. Values greater than or equal to 1 are saturated to just less than 1. A can occasionally evaluate to 
greater than 1 when the control signal is positive. Under these conditions the value is limited to just under 1 . The sine 
interpolator unit 222 is instructed to generate a leading or a lagging sample according to whether A is 1 or 0 respec- 
tively. The spread of sample times covered by the interpolator is one system clock period. 

The rujmerically controlled osdilator 210 operates in a timing loop based on the Gardner algorithm discussed 




so 



7i 



EP0748118A2 



above. Other tirrang recovery algorithms may also be used, such as the Muller and Muller algorithm. A second order 
loop filter 259 is used. A proportional*kTtegral (PI) controller 211 is included in the circuit. Its proportional and integral 
gain constants are selected to give the required damping factor and natural frequency. Preferably a relatively high nat- 
ural frequency is used for initial channel acquisition in order to minimize lock time and insure acquisition. Thereafter the 
coefficients are changed to reduce the loop bandwidth and thereby make It less sensitive to noise and fluctuations. This 
"gear shifting" operation improves the overall system bit erra rate. 

The matched filters 254. 256. prefer abiv square-root raised cosine matched filters, are included because the Gard- 
ner dgorlthm assumes data having no \S\. They carviot be placed before the inta^pdator unit 222 because they have 
hard-wired coefficients designed for T/2 sampled data. As discussed above, tiie Gardner algorithm locks the timing 
sample point using T/2 sanples. The loc^ preferably acquires the sample point such that odd samples are at the zero 
crossing points of the riput data, and the even samples are used as data samples. 

A delay is Imposed between the presentation of system ch)ck rate samples to the sine interpolator unit 222 and the 
appearance of Interpolated samples, according to the following equation 



5 = (system dock period / N); 

N is the nurrt>er of sine interpolation points; 

k :s (integer) interpolation distance, AN; and 

D s constant delay implicit in hardware 
The sine oiterpolator unit 222 is based on a finite impulse response f ater, which is clocked at the system clock rate, 
with the coefficients being selected from a bank of N sets, whereffi each set of coefficient Interpolates a cOfferent delay. 
The interpolation distance output from the numerically controlled oscillator 210 determines whfch bank of coefficients 
are used to generate a given sample, as A varies from 0 to 1 . This can be appreciated with reference to Fig. 5. wherein 
the oval indicators represent interpolation possibilities. Sine interpolatksn is based on the sampling tiieory whk;h shows 
that a signal which has been Nyquist sampled can be reconstructed using sine pulses, equivalent to performing a low 
pass filtering operation in the frequency domain. The output is given by the equation 



Referring to Figs. 15a - 15c. it will be seen that the reconstructed waveform .580 is the sum of all the components, 
representatively shown as lines 582. 584. and 586. There is only one non-zero component at each sample point, as 
shown in Rg. 15b. To make implementation possible, i.e. to make the system causal, the tails of the sine pulse have to 
be trimmed. This introduces ne^igble error. In order to interpolate the value of the signal at a point between known 
samples, it is necessary to sum the contritxttion made at that point by each known sample. The contributions are cal- 
culated based on the amplitude of the samples and the shape of the sine pidse. The FIR filter 250 (Fig. 16) is used to 
calculate and sum the contributtons. The coeffrcients of the filter 250 are calculated based on a system dock rate shic 
pulse. 

As shown in Fig. 16. the finite impulse response (FIR) fitter 250 has a plurality of multipliers 252. each having a 
small read only memory (ROM) 251. The multipliers 252 operate in parallel. Only the ROM 251 for the left-most multi- 
plier 252 is shown for clarity, it being understood that each multiplier is operativeiy associated with a ROM. Many fams 
of memory could alternatively be used. For example in certain applicatkms it may be desirable to operate the receiver 
under control of a mk:roprocessor (not shown) and the memay ROM 251 could be realized as a RAM, with program- 
mable values. Microprocessor interfaces for sanrpling receivers are well known, and wiH not be further discussed 
herein. The coefficients for each delay phase is stored within the ROMs 251 . and an appropriate coefficient is selected 
in the ROM 251 by addressing logic 249 according to tiie interpolatk)n distance provided by the numerically controlled 
oscillator 210. The addressing anangement for the ROMS 251 is conventional. The filter 250 includes a shift register 
257 having pliffality of tapped positfons 258. 

The coeffidents that are stored in the ROMs 251 are based on a system clock rate sine pulse. The zero delay coef- 
fidents have tiie sine pulse centered at zero, and model the latest signal amval when A = 0. The maximum delay coef- 
fidents correspond to a sine pulse centered at (N - 1)/ N system clock periods and are used whenever A > (N - 1)/N . 

Raised cosine sine pulse coefficients, preferably with a 0.35. are generated using the general formula 



delay = D + k5 



where 
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c(x) = [Maths] 
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5 



where 



10 



= = cT + n(I) 



and 

c is the coefficient number (e.g. -2. -1. 0. 1 , 2); 
IS n is the interpolation distance 0...N-1. 

The tap values for a 6 tap. 8 phase interpolator are given in the following table. Row 8 is not implemented. 



Table 1 



20 



25 



30 



6-tap. 8-phase interpolation - FIR filter coefficients 


bank (n) 


A 


Co 












0 


0^A<0.125 


0 


0 


1.000 


0 


0 


0 


1 


0.125^A<0.250 


0.033 


-0.093 


0.974 


0.127 


-0.042 


0.014 


2 


0.250 ^A<0.375 


0.053 


-0.150 


0.895 


0.281 


•0.089 


0.030 


3 


O.375€A<0.50O 


0.061 


-0.171 


0.772 


0.450 


-0.132 


0.046 


4 


0.500 ^A<0.625 


0.057 


-0.163 


0.619 


0.619 


-0.163 


0.057 


5 


0.625SA<0.750 


0.046 


-0.132 


0.450 


0.772 


-0.171 


0-061 


6 


0.750 SA<0.875 


0.030 


-0.089 


0.281 


0,895 


-0.150 


0.053 


7 


0.875^ A< 1.000 


0.014 


-0.042 


0.127 


0.974 


-0.093 


0.033 


8 


not used 


0 


0 


0 


1.000 


0 


.0 



The data generated by the nujttipfiers 252 are summed in an adder unit 253 and output as interpolated data 
The operation of the sine interpolator unit 222 can be further appreciated with reference to Figs. 5. 14a. 14b and 
40 1 5a - c. wherein two exerrplary delays apply. For the examples it is assumed that the interpolator has a 6-tap filter and 
8 phases. In actual practice the number of taps and phases are selected according to the application and the resolution 
desired. In Fig. 14a, the coefficients for Bank 0 is shown. Only the center tap coefficient is non-zero. Thus the output 
data is based solely on the contents of the analog-to-digital value in the position 255. multiplied by a value memorized 
in its associated ROM (not shown). In Fig. 14b a longer delay is involved, and bank 7 of the banks 0 - 8 has been 
45 selected, with the coefficients as shown. For example the coefficient of the center tap has a value of 0.974. Interpolation 
makes use of the known sine pulse shape of the received data. 

Preferably the numerically controlled oscillator, sine interpolator, and ttie ioop filter comprise an integrated semi- 
conducta drcu'it. which may be a CMOS circuit. 

The accuracy of interpolation in the unit may optionally be inaeased by incorporating a level of linear interpolation. 
50 as shown in Rg. 17. which illustrates a linear interpolation unit 267. This performs linear interpdation on sine interpo- 
lated values. As shown in Fig. 5. a required interpolation point 261 is bracketed by preceding and succeeding sine inter- 
polation points 265 and 263 respectively. Linear interpolation to detei'mine the value of point 261 is performed based 
on the points 265. 263. In Fig. 1 7 
8 = (system dock period / N); 
55 N is the number of sine interpolation points; 

k = (integer) Interpolation distance. AN; 
D = constant delay implicit in hardware; and 
f - (fraction) interpolation distance, AN. 
Each system clock sample is sine interpolated at interpolation points K and k + 1 in sine interpolation units 266. 
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268. The interpolation results are multiplied by 1 -f, and f in multipliers 269. 264 respectiveiy. and the result oon*)ined in 

an adder 270. An interpolated sample value is output 
^ Referring to Figs. 1 1 and 18. the matched filters 254. 256 are implemented as finite impulse response filters, and 

are enabled by the valid signal. Strc^e 1 70. that is generated by the interpolator unit 222. An exemplary filter 290 is illus- 
5 trated in Fig. 1 8. Operation of the shift register 280 is enabled with the valid signal 1 70. One output sample is generated 

for each valid input sample. ThfiJElRo o^jcients Co..Cn are c alculated for a T/2 FIR square-rool rajsed cosine filter 
^ assuming an excess bandwidth a = 0.35. T he use of the valid strobe signal 170 emulates docking with a systSii cic^ 

of T/2. although within the filter hardware is actually being clocked at a faster rate, that of the system clock 1 20 (Rg. 1 0). 

^o Carrier Recovery 

The carrier recovery loop Is explained initially with reference to Figs. 1 9 and 20. which illustrates a Costas algorithm 
phase enror estimation section 315, a second order loop filter 320. a numerically controlled oscillator 310. and a digital 
derotation circuit 317. This circuit tracks any frequency errors and phase drift in the external modulation and denxxiu- 
75 lation chains. Preferably an add'rtional adaptive loop circuit 319 operates according to the least-mean-square (LMS) 
algorithm to adaptively estimate demodulation phase ndse en'ors caused by hum and jitter. 

Sin e and cos 0 control the derotation circuit 317. They are generated by use of a look-up table stored in a ROM 
(not shown). The design of trigonometric look-ip tables is well known. 

The derotator rotates the Input data by 6. Given that (l.Q) represents a vector of amplitude (1^ + and argument 
20 tan'^l/Q) ::=<!>. Thus l-sin<|> and Qscos(^. and we require derotated l = r = sin((|) + e) and Q* = cos((|) -i- 6) . 

r s sin^ cose - cos^ sine = I oosO - QsinO; 



25 and 

0* = Qcose + 1 sine. 



30 This is implemented in the network of multipliers and adders shown in the derotation circuit 317. The Costas phase 
enx>r estimation section 315 completes the k)op. 

The derotator 31 7 in cooperation with the phase estimation section 31 5 is also used to correct phase noise and jit- 
ter. This jitter is tracked by an LMS adaptive estimate of the phase error. Referring to Fig. 20, the derotated I and Q val- 
ues, represented as fixed point numbers with a fractional part are sliced in slicers 332. 334 respectively to the nearest 

35 legal constellation value. For QPSK this will be +1 or -1. The difference between tiie derotated value and tiie sliced 
value is dbtained in subtracters 336, 338. and forns the error. The I and Q enror values are converted into a angular 
error estimate e en-or. In the case of QPSK nxxJulation, tiie e error is obtained from a switching network contained in 
angulator 331, according to table 2 below. The output of the angulator 331 is an adapted LMS estimate of the phase 
jitter or hum en-or 9 estimate. The phase error estimation circuit disclosed in our copending plication. Ser. No. 

40 08/481 . 1 07. incorporated herein by reference, can also be used, as may many other phase error estimatbn circuits, for 
example a circuit embodying the Costas algorlthra 

The LMS algorithm and its sign variant is well known, and will not be further explained herein. It is discussed, for 
example, in DigM Communication, Second Edition, by Edward A. Lee and David G. Messerschmitt. Kluwer Academic 
Publishers, Chap. 11. 

45 The adaptive algorithm has been modified slightly from the standard LMS algorithm in that the e estimate has been 
given a leak. Normally leak is zero, but every Nth cycle it is -(sign(e estimate)). This prevents the e error from growing 
beyond operational limits. 

The Costas loop locks with the constellation points on the axis - I.e.. (1. 0), (0, 1). (-1, 0), (0. -1). Thus tine error, in 
the exanple given, as shown in Fig. 24, may be estimated as Oerror = sin*^ (lerror) . which approximates to eerror s 
50 lerror. Similarly for tiie otiier constellation points eerror is either + or - lerror or Qerror as shown in the table below. A 
geometric representation of eerror is shown in Fig. 24. 



55 
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Table 2 



Approximation of e error 


Constellation 


Oeror 


0=1 


1=0 


lerror 


0=0 


!=1 


-Q error 


0=-1 


1=0 


-lenror 


0=0 


l=-1 


Oen-or 



Error calculation for QAM modulation schemes is more complicated. 

15 As in the case of the timing recovery control loop disclosed hereinabove, the proportional and integral gain con- 
stants PI controUer 321 in the second order loop 320 start off with wide bandwidth values to minimize acquisition time, 
and are shifted to a lower bandwidth loop set of values to optimize system bit error rate once lock has been achieved. 
The values selected can be readily selected in accordance v^h the requirements of a particular application. 

After a channel change there may be a significant frequency error which has to be determined before phase can 

20 be acquired. It preferable to implement a frequency-lock-loop (i.e.. one where the error signal is proportfonal to the fre- 
quency error) or a frequency sweeping acquisition scheme to be able to acqure the initial frequency error. A circuit that 
achieves a lock is explained with reference to Rg. 21 . which is associated with the frequency and phase lock loop circuit 
321.. Should the proportional integral loop 320 be unable to lock onto the frequency of the received intermediate fre- 
quency signal, the circuit of Rg. 21 allows the numertcaDy controlled oscillator 310 to "hop" from one frequency to 

25 another at disaete intenals to search for the canier of the incoming signal. The higher-order bits of the output of the PI 
controller 321 , referenced 458 in Fig. 21 . are combined with a hop input 450. taken from a state machine 461, and sub- 
mitted to the hopping adder 414. The adder 414 outputs a frequency offset signal 452 which is accepted by the numer- 
ically controlled oscillator 310. 

30 Second Embodiment 

Many forms of modulation are known to the art which do not generate both the in phase and quadrature compo- 
nents upon demodulation. For example, vestigial sideband (VSB) modulation is achieved by amplitiKle modulating a 
pulsed baseband signal, and suppressing a redundant skieband of the anplitude modulated (AM) signal, in order to 

35 conserve bandwidth. Usually the lower skJeband is suppressed. In the digital form of VSB. a digital pulse amplitude 
modulated (PAM) agnal is employed. The alternate embodiment of the invention disclosed herelnbelow with reference * 
to Rg. 22 is suitable for the reception of VSB signals, as well as many other modulation schemes. As in the first embod- 
iment, the output of an anatog-to-digitaJ converter 560 is applied to a sine interpdata unit 522. which is followed by a 
matched filter 552. and a timing recovery circuit 525. The details of these components are the same as for the first 

40 ennbodiment. and need not be repeated. The derotalion circuit 550. and the earner recovery circuit 555 have the same 
structure as in the first enixxliment. However the derotater circuit 550 requires a quadrature input, whfch must be gen- 
erated, as it is lacking in the sampled demodulated signal that is output by the analog-to-digital converter 560. H is pos- 
sible to operate the timing recovery circuit 525 in accordance with the Gardner algorithm with only the in phase 
component, in which case the error signal given above will be 

45 

error(r) = l[r.|][!(r)+l(r.T)] 



so where 

I is the in phase output; ' 
T is the syntx)! period; and 
r is the sample time of the even sanple. 
The 0 input is developed by a Hilbert filter, shown in Fig. 23. The Hilbert filter has an irrpulse response and a trans- 
55 fer function given by 
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H(j.<D)s=-jsgn(«D) . 



The Hilbert filter is an el even tap Fl R filter, which has been implemented in much the same way as the Fl R filter 290 
5 (Rg. 1 8). The filter is organized as a pluraOty of cells operating in series. accorc£ng to the length of the filter. Or^ cell 
782 is illustrated in Rg. 23. it being understood that the other ceils are structurally identical. In order to reduce hardware, 
the multiplier 786 is shared among the coefficients and the taps *n a data shift re^ster 783. 

The multiplier-accumulator unit 705 of the cell 782 will be described. The data sNft register 783 conprises regsters 
71 1, 712, 713. and 714. and is clocked at T which by way of example is 133 ns. The outputs from the shift registers 
10 710-713 therefore only change only every 133 ns. The cell 782 could be implemented by associating a muJtqslier with 
each of the registers 710-713. for a total of 4 multipliers. However because the multiplier 786 can operate in only 33 ns. 
T/4 the cell has been designed to have one multiplier 786 which is switched by switch 710 between the four data regis- 
ters 71 1-714. Four coefficient registers 720-723 are provided to supply the multiplier 786. Of course it is also required 
that the coefficient registers 720-723 also be switched, indicated by switch 724 In Rg. 23. TTie filter structure requires 
IS that the cell output be formed according to the equation 

CCout=EDnCn 
n»0 

20 

where 

CCout is the cell output: 

On is the contents of the nth data shift register; and 
25 Cn is the contents of the nth coefficient register. 

CCout is accumulated using the adder 726. The individual outputs CCout of the units 705 Is latched, and 

summed in an adder tree 727. As the multiplier requires the largest area of each cell, a large amount of chip area has 
thus been conserved. 

30 Third Embodiment 

A third embodiment of the invention is disclosed herein with reference to Fig. 25. Its construction is similar to the 
first embodiment htowever, having reference to Fig. 10, wftich has been discussed in connection v«th the first embodi- 
ment, It will be noted that the derotater 150 is disposed following the matched f ilters 254. 256. This arrangement has 
35 the advantage of using relatively inexpensive hardware in the derotator. which can be docked at T However the signals 
applied to the matched fitters 254. 256 are affected by constellation rotation and frequency en^ors. and hence the filtered 
output will not perfectly restae the source pulses. In Rg. 25 the derotater 652 is disposed intermediate the interpolator 
622, and the matched filters 754, 756. The derotater 652 is now reqiired to be clocked at T/2; however the signal pro- 
duced by the matched filters 754, 756 is a more accurate restoration. 
'40 While this invention has been explained with reference to the stoxjcture disclosed herein, it is not confined to the 
details set forth and this application is intended to cover any mo(§f icatfons and changes as may oome within the scope 
of the folfowing dairns: 

Clalnis 

45 

1 . A signal processing apparatus for processing signals that are sampled by a sampler operative at a sampling rate, 
the signals having a received symbol rate, the apparatus comprising: 

a clock, operative at said sarrpling rate; 
50 a first numerically controlled osdilator operative at periods T that are initially equal to a nominal baud rate 

divkied by said sampling rate; 

a sine interpolator receiving samples at said sanpling rate; and 

a loop filter, coupled to said sine interpolator and said first numerically controlled oscillator and having an out- 
put responsive to a difference between said periods T and said received symbol rate of sakl sampled signals; 
55 wherein said first numerically controlled oscillator is responsive to said loop filter and generates an output sig- 

nal that is representative of an interpolation distance between succeeding samples, and sakl interpolator inter- 
polates said received samples according to said interpolation distance, and produces an output signal 
representative of said interpolated samples. 
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2. The apparatus according to claim 1 . wherein said sine interpolator performs a plurality of sine interpolations that 
precede and follow a required sine interpolation point, further comprising a linear interpolator that performs linear 
interpolation on said plurality of sine interpolations. 

5 3. A circuit for processing nxxiulated signals, comprising a semiconductor integrated carrier recovery circuit operative 
to control a demodulator, the circuit comprising: 

a second numerically controlled oscillator; 

a digital derotation ctrcuH responsive to said second numerically controlled oscillator and accepting an in phase 
10 component and a quadrature component of sampled signals: 

a phase error estimation circuit, coupled to an output of said derotation drcuH; and 
a loop filter coupled to an output of said phase error estimation circuit; wherein said second nimierically con- 
trolled oscillator is responsive to said loop filter. 

15 4. The circuit according to daim 3. further comprising a drcuit for adaptively estimating the phase error according to 
a least-mean-square algorithm, the circuit including: first and secortl slicers. accepting a derotated in-phase value 
andaderotated 

quadrature value respectively: 

first and second subtracters, for respectively determining first and second differences between said dero- 
20 tated in phase value and said sliced in phase value, and between said derotated quadrature value and said sliced 
quadrature value: and 

an angulator. accepting said first and second differences and outputthg a phase enror estimate. 

5. A signal processing apparatus for processing modulated isignals that have a modulation carrier frequency, tiie 
25 apparatus oomprisoig: 

a demodulator; 

a sampler operative at a sampling rate on an output of said demodulator; and 
a carrier recovery circuit operative to control said demodulator in accordance with 

30 

said modulation carrier frequency* the drcuit comprising: a second numerically controlled oscillator; 
a digital derotation circuit responsive to said second numerically controlled oscillator and accepting an in 
phase conponent and a quadrature component of sampled signals: 
a phase enror estimation drcuit, coupled to an output of said derotation drcuit: and 
35 a loop filter coupled to an output of said phase error estimation circuit: 

wherein said second numerically controlled osdilator fs responsive to said loop filter: 

wherein said sarrpler and said carrier recovery circuit are integrated in a semiconductor integrated drcuit. 

40 

6. A signal processing apparatus for processing modulated signals that are sampled by a sampler operative at a sam- 
pling rate, the signals having a received symbol rate, the apparatus comprising: 

a clock, operative at said sampling rate; 
45 a first numerically controlled oscillator operative at periods T that are initially equal to a nominal baud rate 

divided by said sampling rate; 

a sine interpolator receiving samples at said sampling rate; 

a first loop filter, coupled to said sine interpolator and said first numerically controlled osdilator having an out- 
put responsive to a difference between said periods T and a transmitted symbol rate of said sampled signals; 
50 wherein said first numerically controlled oscillator is responsive to said first loop filter and generates an output 

signal that is representative of an interpolation distance between succeeding sanples. and said sine interpo- 
lator interpolates said received sanples according to said interpolation distance, and produces an output sig- 
nal representative of said interpolated samples: and 
a carrier recovery circuit, comprising: 

55 

a second numerically controlled oscillator; 

a digital derotation drcuit responsive to said second numerically controlled oscillator and accepting an in 

phase component and a quadrature component of sampled signals: 

a phase error estimation drcuit. coupled to an output of said derotation drcuit; and 
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a \oop filter coqpled to an output of said phase error estimation drcuit: wherein said second numerically 
controtled oscillator is responsive to said loop fater; 

wherein said first and second numerically controlled oscillators, said sine interpolator, said first and second loop fil- 
ters, and said dgital derotation circuit are integrated in a semiconductor integrated drcuit. 
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